SO(N) symmetries in the two— chain model of correlated fermions 
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The two-chain model of correlated fermions, appropriate e.g. for doped spin ladders or carbon 
nanotubes, is investigated in the weak-interaction limit. It is shown that the model scales to 
situations with high internal symmetries, generically S0(6), whereas on critical lines there is SO(5) or 
SO (3) symmetry. Consequences for the level spectroscopy of the model are discussed. In particular, 
the "7r-mode" is shown to be sixfold degenerate and contains both spin and charge excitations. 



Models of two coupled chains of interacting fermions 
are of interest for a variety of reasons: they are directly 
applicable to interesting experimental systems like doped 
spin ladders |jj or carbon nanotubes ||, and they also 
provide a first step for the understanding of the transition 
from the well-understood single chain case to the area 
of two-dimensional correlated fermions. || In particular, 
both analytical and numerical work show that 
the two-chain system with purely repulsive interactions 
is a d-type superconductor whereas the single chain is 
antifcrromagnetic. This provides a potential explanation 
for the observation of superconductivity in a doped spin 
ladder, Q and is probably the best-established case of 
correlation-induced superconductivity. 

In recent work fl(J it has been shown that the Mott- 
insulating half-filled two-chain model has a very big 
SO(8) symmetry in the weak-correlation limit, leading 
to a number of interesting predictions. Away from half- 
filling the SO (8) symmetry is broken and the system 
becomes conducting. Here I will study the symmetries 
of the conducting state away from half-filling and show 
that generically a still large SO (6) (=SU(4)) symmetry 
remains. On critical lines the symmetry is lower, SO (3) 
or SO (5). There are exact solutions corresponding to 
these symmetric models, jll],[l2| and from those it is pos- 
sible to obtain quite detailed information about the ex- 
citation spectrum and in particular the appearance and 
location of bound states. The present work also provides 
a microscopic extension of the SO (5) model of correlated 
fermions |l3| and the associated ladder models. jL4j 

I consider the standard two-chain model of interacting 
fermions. [|]-[| After linearizing the electron dispersion 
around the Fermi energy, the kinetic energy density takes 
the form , 



h = -iV F i^Rsg^R 
s,q—0 : 7T 
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Here q = (it) for the (anti-)bonding bands split by the 
interchain hopping integral tj_, s is the spin, R (L) refers 
to the right- (left-)going fermions, and tpR(L)sq are the 
corresponding fermion field operators. 

Important interaction effects arise from processes 



where both the two incoming and the two outgoing parti- 
cles are all at the Fermi energy. The Hamiltonian density 
of the interaction then is 



hint = 9l,abcd^ Rsa lpl tb lp Rtc i> 
+ 92,abcd^ Rsa ^ Ltb ^ Ltc ^ 



Lsd 



Rsd 



(2) 



These terms describe backward (51) and forward ((72) 
scattering between particles at the Fermi surface. Sum- 
mation over the spin indices s,t and the band indices 
a, b, c, d is understood. The interactions allowed by mo- 
mentum and energy conservation are 51,0000 = 9\,^-kw = 

911 , gifiOirir = ffl,7nr00 = .912, 9l,0ir0ir — ffl^OwQ = <?13, 
52,0000 = 52,7rir7TTr = 521, 32,0irir0 = ,92,7r007r = 522, 
52,007T7r = <72,7T7rOO = <?23- 

The one-loop renormalization group equations for the 
coupling constants under a change of energy cutoff A — > 
Ae~ l have been found previously Ala] and are given by: 
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Here g> ? = 2irv F ^ ^ 

Abelian bosonization is used to write the model in 
terms of symmetric and antisymmetric charge (y = p±) 
and spin (v = cr±) fields <j> v and their conjugate mo- 
mentum densities n„. The total Hamiltonian becomes 
the sum of a term bilinear in boson fields (i?&) and a 
term containing nonlinear interactions H n \, with Hb = 



M v ) ILi + -= (1 - M u ) (dxfa)' 



and M a ± = (511 ± giz)/2-KVF ,M P ± = M a ± - (521 ± 
522) /ttuF' The interaction term is 



dx{(g 2 3 ~ 512) cos 28 p _ cos 29 a . 
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Another transformation, which will be more generally 
useful is obtained by refermionizing the bosonized Hamil- 
tonian, eq.(Q), by introducing three fermion fields ^1.2,3 
corresponding to the (ft—), and (c— ) boson fields. 

After a duality transformation </> p _ «-* and writ- 
ing the Dirac fermions in terms of hermitian (Majorana) 
fields x as ijj a i = Xa,2i-i + *Xa,2« the Hamiltonian in this 
subspace then takes the form 



H = H + 2g / dx 



SN=f> \ 
XRiXLi 



(7) 



FIG. 1. "Phase diagram" in the g\ — gi plane 
(gu = gi,gn = 52), indicating the dominant fluctuation. 
CDW: charge-density wave, SCd: d-wave superconductor, 
OAF: orbital antiferromagnet, SCs: s-wave superconductor. 
On the critical lines g-z = and g\ = 2g2 additional critical 
fluctuations occur. 

+ 523 cos 29 p - cos 2<p a - + gu cos 20 CT+ cos 2<^> CT _ 
+ 312 cos 20 p ^ cos 2<j) a+ - 513 cos 2</> CT+ cos 20 a J$) 

where 9 V is the field dual to 4> v : U u = d x 9 u /ir. 

It is clearly interesting to understand the energetics 
and quantum numbers of the low-lying excited states 
of the model, in particular in order to understand the 
relation to numerical or experimental results. I start by 
investigating the generic regions of the phase diagram, 
fig^, i.e. away from the critical lines g\ = 0, g\ — 2g 2 . In 
particular, in the SCd region, important for Hubbard- 
like bare interactions, a numerical investigation shows 
that the asymptotic behavior near the divergence (I — > l c ) 
is given by g%j{l) = gfj/Hc — I), with the values 

ff 2 °3 = l/4, 9° 13 = 0, 
3/16, g° 22 = 1/16 . (5) 

One can note that eq. is an exact solution of the renor- 
malization group equations (^) for a special set of initial 
conditions. For more general interactions, in particular 
the Hubbard model, the solution of the renormalization 
group equations tends asymptotically for I — > l c towards 
this solution, i.e. eq.(Q) is the "fixed ray" towards which 
all solutions of the renormalization group equations scale 
in the SCd part of the phase diagram. 

Interestingly, for the parameter values eq.(|5|) the model 
exhibits a high but somewhat hidden symmetry: after 
the particle-hole transformation on the left-going branch 
^Lsq i^L-sq takes the C7(4) invariant form 
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H = H -g J dx{^p R p L + ^ Rsq i> Lsq i>l 



Ltp4>Rtp)>( 6 ) 

where p a = "4>tsq4>asq- This is the chiral Gross-Neveu 
model, for which an exact Bethe ansatz solution exists. 

Pi 



This is the well-known Hamiltonian of the Gross-Neveu 
model [|l(| which has a SO (6) symmetry. Due to the 
equality of the Lie algebras so(6) and su(4), the equiv- 
alence of the chiral SU(4) and the (non-chiral) SO(6) 
models may seem quite natural. One should however no- 
tice that the nonchiral model has a double-degenerate 
symmetry-broken ground state, |16 1^] whereas there is 
no symmetry breaking in the chiral model. This will lead 
to interesting consequences for the quantum numbers of 
the excited states. It should be emphasized that eq.(Q) 
describes the massive part of the spectrum and that the 
total charge mode p+ remains massless in all cases and 
therefore contributes to correlation functions in an im- 
portant way. 

We can now use known results for the excitation spec- 
trum of the SO(N) models ]Tl| , p^[ to obtain results for 
the excitations of the two-chain model. For the generic 
SU(4) case the elementary excitations are single fermions 
and holes with mass gap mo oc e^ 1 ^ 9 . Due to the 
particle-hole transformation leading to eq.(Q), right- 
going fermions transform as the fundamental 4 repre- 
sentation, whereas left-going fermions transform as the 
conjugate 4. In the SO (6) language these fermionic ex- 
citations are kinks separating the two degenerate ground 
states, (T7) whereas the SO(6) fermions correspond to 
two-particle excitations in the original model. 

From the single-fermion gap tuq one expects a gap 2too 
in two-particle excitations. This is in fact correct for all 
excitations created by operators of the form i^Rsq^Ltq 1 • 
Consequently, all pairing correlation functions, irrespec- 
tive of their spin (singlet or triplet) or orbital structure, 
decay exponentially, with the same correlation length 
£ = vp/2mo. The one important exception of course is 
the "d-wave" singlet channel, which has purely algebraic 
decay. In the group theoretically language this corre- 
sponds to the decomposition 4 <g> 4 = 10 15, the 1 being 
the d-wave superconductor. 

More structure arises in the particle-hole channel, cre- 
ated by operators of the type V , j? S9 '0Li g " Here the rele- 
vant decomposition is 4 ® 4 = 6 © 10, the 6 being the 
antisymmetric rank-2 tensor. Remarkably, the exact so- 
lution shows that particle-hole excitations in this channel 
are bound and appear at energy \^2mo, rather than 2mo- 
The operators creating these bound states are the spin 
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triplet V'fl s o( (T a)st' ( /'Lt7r + (0 7r ) w hich represents an 2kp 
spin density wave which is out of phase between the two 
chains, and the "orbital triplet" ip]i sq {o'a)qq'ip[ jS qi which 
contains a charge density wave which is in phase between 
the two chains and the orbital antiferromagnet (OAF, or 
staggered flux phase). From the SU(4) symmetry we can 
conclude that the gaps in these six channels are identi- 
cal. It is in particular remarkable that the lowest spin 
and charge gap appear, due to bound state formation, at 
only y/2 times the single-particle gap. One should notice 
that contrary to the triply-degenerate antiferromagnetic 
7r-mode of the SO (5) model of correlated fermions p3[ , 
we here have an extended ir-mode with sixfold degeneracy 
and which also contains charge excitations. The remain- 
ing particle-hole excitations (the 10) all have excitation 
energy 2m . 

I now turn to the critical lines in fig|l[ Both are of Ising 
type. For g\ = the asymptotics of the flow is now given 

by 9u = 0,9°i = -3°2 = 323/ 2 = V 4 - 0nc can now fol- 
low the same logic as before. However, here the self-dual 
combination cos 2<j) a - +cos 29 a - — > XR3XL3 appears, and 
consequently Xa4 decouples and remains massless, to- 
gether with the two boson fields <f> u +- The fermionic 
model for the massive degrees of freedom then is given by 
eq.(0), however now with N = 3. In this case the model 
has an SU(2) fl x SU(2)l x SU(2) symmetry, where the 
first two subgroups represent spin rotations for the right- 
and left-moving electrons (these are independent because 
gij = 0), whereas the last SU(2) acts in an "orbital" 
space spanned by the spinors — (tpR S o, *Prstt) T and 

There is a gap mo in the single-particle (kink) excita- 
tion spectrum. The two-particle and particle-hole spec- 
tra can be analyzed rewriting the p— and <r— fields in 
terms of Ising model variables, jlq] This then shows that 
excitations of orbital singlet type are massless. In the 
particle-particle channel this leads to powerlaw decay of 
d-wave pairing correlations, both in the spin singlet and 
triplet channels. In the particle-hole channel, the pow- 
erlaw decay occurs for CDW W and SDW T correlations 
which are out of phase between the two chains. The ex- 
ponent of the powerlaw is entirely determined by the total 
charge degrees of freedom (p+). For N < 4 the Gross- 
Neveu model does not show two-kink bound states. Jll]] 
Consequently, one easily checks that all other particle- 
hole or two-particle excitations have the same gap 2too, 
and the corresponding correlation functions decay expo- 
nentially with distance. These states are all triplets un- 
der the orbital SU(2). 

Finally, on the second critical line g\ — 2g2 there is 
again a self-duality in the a— sector, but now the o~+ 
degrees of freedom also have a gap. Here g^ = g° 3 = 

.921 = -323 = -l/6,3?2 = 1/3,322 = 0. We again 
find a Gross-Neveu model, but now with N = 5. The 
kink excitations now form a quartet, leading again to a 
gap mo in the single particle Green functions. Arrang- 



ing these quartets as = (tpR^o, V'fij.o, iipRU, i4>Rin) T , 
and *l = (#LTir, iipLiK, ipLio, ipLioV, these transform 
as SO (5) spinors, and the Hamiltonian with these cou- 
pling constants given by then is explicitly SO (5) invari- 
ant. 

Contrary to the SO(3) case, there are now two-kink 
bound states at energy ^/3mo, in addition to unbound 
two-kink states at 2mo. These bound states occur in 
the vector (5) part of the SO (5) decomposition 4 £g> 4 = 
1®5®10. In particular, in the particle-hole channel, the 
scalar (1) excitations are massless and occur in the OAF 
channel, and the bound vector excitations are the same 
as those in the particle-hole channel in the SU(4) case, 
with OAF removed. All other particle-hole excitations, 
occurring in the antisymmetric tensor (10) representa- 
tion, have mass 2mo- In the particle-particle channel the 
same decomposition applies, however now the massless 
scalar corresponds to d-wave singlet pairing, the bound 
vector consists of s-wave triplet pairing and two singlet 
paired states with transverse momentum tt, and the 10 
contains all other pairing states. 

The energy levels in the interesting region < g± < 
232, corresponding to repulsive interactions, and on the 
adjacent critical lines are shown in fig.§ In the SO(N) 
language, these are the discrete energy levels at q = or 
q = 2k p, however in our original model the presence of 
the massless p+ mode transforms these into lower edges 
of continua. Near this edge, the correlation functions 
have singularities of the form (lo — ljq)~ u , where in weak 
coupling v « 1/2 and for the lowest (massless) levels 
wo = and v w 3/2. 

The full symmetries are only exact for asymptotically 
weak coupling. At any finite coupling, corrections to the 
asymptotics of eqs.(^|) will lead to splittings in some mul- 
tiplets. This is schematically indicated for the bound 
multiplets in fig.§. A more dramatic consequence of the 
corrections occurs in the close vicinity of the critical lines: 
for example, in order to obtain the massless CDW^ cor- 
relations of the SO(3) case, the correction terms have to 
form a bound state out of the continuum of the SO (6) 
spectrum and to "pull it down" all the way to mass- 
lessness as the SO (3) critical line is approached (as in- 
dicated by the broken lines in fig.||). This is in par- 
ticular of interest for the case of very weakly coupled 
chains where renormalizations from the one-dimensional 
region E > t± lead to initial conditions for eqs.(^) obey- 
ing gu <gC <72i- A detailed discussion is beyond the scope 
of the present paper. 

Other parts of the phase diagram, fig.|l|, can be dis- 
cussed along similar lines. Away from the critical lines 
one always finds an SO (6) symmetry. Of some interest is 
the SCs region, which contains the Hubbard model with 
attractive interaction which is numerically quite acces- 
sible and might serve as a testing ground. Here again 
a sixfold degenerate low-lying state exists, however now 
composed of the three CDWo.^ states and the spin ne- 
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FIG. 2. Energy levels of the SO(N) models for different particle-hole and particle-particle states in the region < g\ < 2<?2 
of the phase diagram, fig|l|. Horizontal lines indicate discrete energy levels, shaded areas two-particle continua. In the SO(5) 
case, the bound state at y/3rno also appears in the two-particle channel, as discussed in the text. Labels indicate in which 
correlation function the corresponding level appears. The lowest levels in each column indicate massless correlations. TSCd: 
triplet d-wave superconductor; (S)CDW 9 : (spin) charge density wave in phase (q = 0) and out of phase (q = 7r) between 
chains. 



matic (the spin triplet analogue of the orbital antiferro- 
magnet). 

In conclusion, I have shown that in the weak- 
interaction limit the generic doped spin ladder model ex- 
hibits an SO (6) symmetry in its spectrum, weakly broken 
by corrections to scaling. The most direct consequence 
is that there is a sixfold degenerate extended n-mode in 
the d-wave superconducting state which contains both 
spin and charge excitations. This is to be compared with 
the SO(5) model [^3) (where the symmetry is broken to 
SO (3) by doping) and only a triply degenerate spin 7T- 
mode is predicted. Approximately equal charge and spin 
excitation energies have been observed in numerical cal- 
culations |l9| on t-J ladders, providing an indication that 
the present results may be relevant even for strongly cor- 
related systems. To which extent the present results can 
be generalized to genuinely two-dimensional situations 
remains to be studied. 
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